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Abstract 

We show that the alternative discrete Painleve I equation (alt-dPi) has a unique 
solution which remains positive for all n > 0. Furthermore, we identify this positive 
solution in terms of a special solution of the second Painleve equation (Pii) involving 
the Airy function Ai(t). The special-function solutions of Pn involving only the Airy 
function Ai(t) therefore have the property that they remain positive for all n > 0 and 
all t > 0, which is a new characterization of these special solutions of Pn and alt-dPi. 


1 Introduction 

We will investigate the system of two nonlinear equations 

Qn ®n+l V b (1.1a) 

an{hn + hn-i) = n, (1.1b) 
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which is known as the alternative discrete Painleve I equation (alt-dPi) [7]. These equations 
arise, for example, when one wants to find the recurrence coefficients of orthogonal polyno¬ 
mials with an exponential cubic weight; see, e.g. mm- They also arise when one wants to 
compute the recurrence coefficients of multiple orthogonal polynomials with an exponential 
cubic weight [6]. The equations also give relations between solutions of the equation 


^ = V - 2*9 - 2a, 


( 1 . 2 ) 


for different values of the parameter a mm- Equation ( |1.2[ ) is equivalent to the second 
Painleve equation (Pn) 


d^w 

dz^ 


= 210 "^ -j- zw + a, 


(1.3) 


through the scaling w{z) = —2 ^^^y{t), with z = —2^/^t. We note that if we solve ( l.lb| ) 
for a„ and substitute this in (1.1a), then we obtain 


n -|-1 n ,9 

7- 7 -+ 7- 7 - = K-t. 

Ofi T ^n+l T ^n—1 


(1.4) 


The main result in this paper is the following. 


Theorem 1.1. For nonnegative values oft, there exists a unique solution of (1.1) with 
ao(t) = 0 for which On+iit) > 0 and bn(t) > 0 for all n > 0, corresponding to the initial value 

bait) = -A\{t)/A\{t), (1.5) 


where Ai(t) is the Airy function. 


The corresponding result for the nonlinear recurrence (1.4) is 


Theorem 1.2. For nonnegative values oft, there exists a unique solution of (1.4) for which 


bn{t) > 0 for all n>0, corresponding to the initial value (1.5) 


We observe that the value of bi(t) is fixed once ao{t) and bo{t) are known. Specifically 
one has bi{t) = —bo{t) + l/[bl{t) — f]. Hence the solutions of ( |1.4 ) only depend on one free 
parameter bo{t) even though it is a recurrence relation of second order. This is reflected in 
Theorem |1.1| by the initial value ao{t) = 0. It is important to note that the initial value ( 1.5[ ) 
involves only Ai(t), rather than a linear combination of the Airy functions Ai{t) and Bi(t). 

Uniqueness of positive solutions for non-linear recurrence relations was already studied 
earlier for the equation 

Xn{Xn+l + Xn + Xn-l) = ri, (1.6) 

which is known as a discrete Painleve I equation mm and which appears naturally when one 
wants to find the recurrence coefficients of the orthogonal polynomials with the exponential 
weight exp(—a:^) on the real line, the so-called Freud weight. Lew and Quarles [I3], Neva! 


T6] . and Bonan and Neva! [3] showed that there exists a unique solution of ( |1.6 ) with xq = 0 
for which > 0 for all n > 1 and this solution corresponds to Xi = 2P(|)/P(1). This result 


was recently also proved for a family of non-linear recurrence relations generalizing (1.6), see 

m- 
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The fact that there is a unique solution of (1.6), or equivalently for (1.4), means that it 


is not wise to compute this positive solution starting from the initial values ao(t) = 0 and 
bo{t) = —Ai'(t)/Ai(t), because a small error in bo means that one will not be generating 
the positive solution and after some time there will be negative bn{t). In (|^we will give a 
method to compute the positive solution (6n)n>o by means of a hxed point algorithm using 
a contraction acting on inhnite sequences. That method is stable and is very suitable to 
compute the positive solution numerically. 

One can easily obtain the asymptotic behaviour of the positive solution of (1.6). This 
was already proved in [6] , but we repeat the proof here for completeness. 


Corollary 1.3. For the positive solution of (1.6) one has hm„_,.oo = 

lim„^oo bn/n^^^ = 1. 


and 


Proof. The positivity of the a„ and (|l.la) imply that On <b‘^. The positivity of the bn and 


(1.1b) imply < n. Together this implies a„ < r?jaf^ so that a„ < rpl"^. From (1.1a) it 


then follows that b\ < rifl'^ + (u + 1 )^A_ We can then dehne 

Bi = liminf B 2 = limsup 


If we take a subsequence such that bnlnfl"^ —)■ i?i, then from (1.4) one hnds 

B\>2I{B^FB2). 


In a similar way, by taking a subsequence such that bn/n^^^ —)■ h follows from (1.4) that 


B2 — 2 /{Bi + B 2 ). 


Together this implies that i ?2 < B\, and since one always has Bi < B 2 , it follows that 
Bi = B 2 so that bn/n^^^ converges. From (1.1b) it then also follows that anlrpl'^ converges. 
Taking limits in (1.1a) and (1.1b) then easily gives the values of these limits. □ 


In 


This paper is organised as follows 
behaviour in (|3| In (0we discuss the relationship with solutions of Pn (1.3) 


we discuss the uniqueness of and discuss its 

It is well known 


that Pn (1.3) has special-function solutions expressible in terms Airy functions |H HOl [T7] . 


Using these we give explicit expressions for the solutions of (|1.1|) and (|1.4|) in terms of the 
Airy function Ai(t) in 


The special-function solutions of Pn involving only the Airy 
function Ai(t) therefore have the property that they remain positive for all n > 0 and all 
t > 0, which is a new characterization of these special solutions of Pn and alt-dPi. 


2 Proof of the uniqueness 


In this section we will prove Theorem 1.2 A positive solution {bn)n>o implies that oq = 0 
and On > 0 for n > 1 by using (1.1b), so that Theorem O follows as well. The idea of the 
proof is to construct a mapping T on the space of positive sequences (xn)n>o and to show 
that this is a contraction on a complete set in this space. The unique hxed point will be the 


desired positive solution of (1.4). 
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Let {xn)n>o be an infinite seqnence with > 0 for all n > 0. We define a new seqnence 
((Ta:)„)„>o implicitly by 


n + 1 


n 


and 


(T^)n “1“ 2^71+1 (T^)n “I" Xji—1 

1 


= {Tx)l-t, n>l 


(Tx)o + xi 
Observe that (Ta;)„ is a solntion of 

n + 1 n 


= [{Tx)or-t. 


y + Xn+I y + Xn-l 


= y -t. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


If Xn+i > 0 and Xn-i > 0, then the left hand side of (2.3) is a positive and decreasing 
fnnction of y E (0, cx)) and the right hand side of (2.3) is obvionsly increasing on [0, cx)) so 
that (2.3) has a nniqne positive solntion and therefore for n > 0 we set (Ta:)„ to be the 
positive solntion of (2.3). Note that the left hand side reqnires Xn+i and Xn-i for n > 1 bnt 
only reqnires Xi for n = 0. Eqnation (2.3) corresponds to a qnartic equation in ?/ for n > 1 
and a cubic equation if n = 0. There is a negative solution between —Xn+i and —Xn-i when 
n > 1 and two complex conjugate solutions or two real negative solutions which will not be 
used in the remainder of this paper. 

We will now give some properties of this mapping T. 


Lemma 2.1. If Xn > 0 for all n > 0, then 0 < (Ta;)„ < Bn{t) with 
Bn{t) = {n+ \ 


t^ 



1/3 


1)2 

21 

+ 

1 + ^ 


t^ 


-I 1/3 


27{n + \Y 


(2.4) 


Proof. Since Xn+i > 0 and Xn-i > 0, it follows from (2.1) 


Xn = (Ta;)„ = \ t + 


n + 1 


+ 


n 


(Tx)„ + Xn+l (Ta;)n + Xn-l 


< \ t + 


2n + l 
{Tx)n 


and therefore 


(T(r)„| [(T(r)„] ^ < 2n + 1. 

The roots of the equation y{y‘^ — t) = 2n + 1 can be written as: 



-1 1/3 


+ 

1 + ^ 




1 1/3' 


y2± = 


P 

27(n + 4)2 


1/3 


1 + 4 / 1 — 


P 


27(n+ i)2 


11/3 


27(n + 1)2 
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(2.5) 


In fact Vi = {n+ ^y^^Q{w), with w = ^t{n + \) where 

Q{w) = — a/ 1 — j + a/1 — 

Notice that for to > 1 it follows that Q{w) = 2^/wcos [| arccos(tc“^/^)]. Hence, Q{w) is 
always positive for any real value of w, in fact it is the unique positive solution of the equation 

Q{w) \Q‘^{w) — 3tc] = 2. (2.6) 

As a consequence, yi is positive for any real value of t (and any integer n). The roots y 2 ± are 
either a pair of complex conjugate numbers, when (2n + l)^ > or two negative numbers, 
when (2n + 1)^ < The case where eventually (2n + 1)^ = would produce a positive 
single root yi and a double negative root yi-. Hence, no matter the value of t, there will be 
a unique positive root, which is yi = Bn(t). Now the positivity of (Tx)„ yields the required 
bounds. □ 


Let us write Bn{t) = i?(n+ |,t) with R{z,t) = (t/(3^^/^)), where Q{w) is the 

function given by (2.5|). 


Lemma 2.2. The function 
R{z, t) = z^^^ 


R 


1 1/3 


+ 


1 + f 1 


R 

TTR 


1 1 / 3 ' 


is increasing and concave in z on (0, oo) and increasing in t on M. 


Proof. Since the function Q{w) is the unique positive root of the equation (2.6) it follows 
that Q'^iw) > 3w and therefore Q‘^{w) > w for any tc > 0. Now, by differentiating (2.6) 
once we obtain 

Q'{w) [Q^{w) -w] = Q{w) > 0, 
and consequently Q'{w) > 0 and also 


Q{w) — 2wQ'{w) = 


Q{w) 


Q‘^{w) — w 


[Q‘^{w) — 3w] > 0. 


A second differentiation of (2.6) gives 


Q"{w) [Q^{w) -w] = 2[l - Q{w)Q'{w)]Q'{w), 


which can be written as 

Q”{w) [Q‘^{w) -w]^ = -2wQ{w), 

and therefore we conclude that Q"{w) < 0 for any tc > 0, and Q"{w) > 0 if tc < 0. 
Now, from the dehnition of the function R{z,f), observe that 


dR 

dz 


1 

3^ 


Q{w) — 2wQ'(w) 




1 Q{w)[Q‘^{w)-3w\ 

3z^R Q‘^{w)—w 
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and also that 


d^R _ 2 

dz'^ 9 z^R 

2Q"{w 


= < 9wz^R 


w=t/{3z2/3) 
_24/3/(9^5/3)^ 

^ , OR d^R , OR 

Consequently, > 0, —^ < 0 and —— > 0. 

oz oz^ at 


Q{w) + 3wQ'{w) + 2w^Q"{w) 

if t > 0, 
if f = 0. 


w=t/(3z^/^) 


□ 


Lemma 2.3. Ifxn < Bnit) for alln > 0, then (Ta;)„ > CiBn{t) for alln> 1 and (Ta;)o > Cq, 
where Cq = 0.68554389 and Ci = 0.6379714. 


Proof. Equation (2.1) easily gives 




n + 1 


+ 


n 


{i:x)n + Bn+iit) iTx)n + Bn_i{t) 


for n > 0. 


(2.7) 


With Lemma 2.2, we readily observe that for z e [0, oo) the function f{z) = l/[y + R{z, t)] 
is such that 

df 1 dR 

< 0 , 


whilst 


d^f 


dz [y + R{z,t)]‘^ dz 

1 fdPy _ 1 


d‘^R 


> 0 . 


dz'^ [y + R{z,t)]^ \dz J [y + R{z,t)]'^ dz"^ 

Hence f{z) is convex on (0, cxd), so that 

A/(^i) + (1 — A)/(z 2) > /(Azi + (1 — A)z2), a G [0,1], 

whenever 0 < < ^2 < oo. For n > 1 we choose Zi = n + |, ^2 = | and A = ^, then 

this gives for y = (Ta;)n 


> 


(Ta;)n + i? (n + |,f) {fTx)n +R{n - \r) {Tx)n + R {n + 
which easily leads to 

n + 1 n _ 2n 


for n > 1, 


> 


(Ta;)n + i? (n + |,f) {fTx)n +R{n - \r) {Tx)n + R {n + 
Combined with ( |2.7| ), and also because Bn{t) = R{n + this gives 

2n 


for n > 1. 


(Ta;)„ - f > 


(Tx)„ + Bn{t) 


, for n > 1. 


and therefore 

[iTx)l - t] [(Ta;)„ + > 2n, for n > 1. 
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With c = (Tx)„/-B„(t), the latter inequality 


(c + 1) c^{Bn{t)) —tBn{t)>2n, for n > 1, 
which can be expressed as 

(c + 1) (c^ — l)B^{t) + {2n + 1) > 2n, for n > 1, 


because tBn{t) = B^{t) — {2n + 1). Since Bn{t) = (n + Q{wn{t)), for n > 0, with 

Wn{t) = t/[3{n + it follows 


(c+ 1) (c^ - 1) Q{Wn{t)) 


1 3 


+ 2 ^ > 


An 


2n + 1’ 


for any n > 1. 


( 2 . 8 ) 


For f > 0, the sequence [wn{t) = t/[3{n + is monotonically decreasing so that 0 < 

Wn(t) < wi{t), n > 1. Since Q(z) is a positive, increasing function (with hni 2 _,._oo Q{z) = 0, 
lim^^+oo Q{z) = +oo), it follows that 


Q(0) < Q{wn{t)) < Q{wi{t)), 
for any f > 0 and n > 1, where Q(0) = 2}!'^ and 

1/3 


1/3 


QK W) = - \/i - + 1^1 +-^1-^^(f) 

Since 0 < (Tx)„ < B^it) implies 0 < c < 1 and therefore (c + l)(c^ — 1) < 0, then 


2c2(c+l) > (c+l){(c'-l)[Q(w„(f))]^ + 2} > 


An A 

- > 

2?7, + 1 3 


Hence, c is such that c^(c+ 1) > The function 2c^(c+ 1) is monotone increasing on [0,1], 
and this implies c > ci, where ci is the positive root of c^(c +1) = |, which is approximately 
Cl = 0.6379714. 

For n = 0 we use (2.7) to obtain 

[(Ta:)o]2-t>l/[(Tx)o + Hi(t)], 


so that [(Tx)o + Bi{t)] {[(Tx)o]^ — t} > 1. With d = (Tx)o/i?i(t), then we have 

(d + 1) [d^Bfit) - tB^it)] > 1. 

With the identity tR (|,t) = Bf{t) — 3, the latter inequality becomes 

(d+1) [{d^-l)Bf{t) + 3] > 1. 

which, because of the fact that Bi{t) > 3^^^ and 0 < d < 1, implies 

(d + 1) [(d^ - 1)3 + 3] > (d + 1) [(d^ - l)B^{t) + 3] > 1. 

The function d^(d + 1) is monotone increasing on [0,1], and this implies d > cq, where 
Co is the positive root of d^(d + 1) = which is approximately Cq = 0.47533, so that 
(Tx)o > 0.47533 Hi(t) > 0.47533 x 3^/3 = 0.685544. □ 
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If we use the norm ||a:|| = sup„>Q with as given in (2.4), then Lemma 


2.1 implies that ||Tx|| < 1 whenever {xn)n>o is a positive sequence. Lemma 2.3 implies that 
||T^x|| > Cl. We are interested in the iterations T^x for k —)■ oo, so without loss of generality 
we may assume that ciBn{t) < x„ < Bn{t). The mapping T is then a mapping from positive 
sequences in the unit ball Bi = {(x„)„>o : ||x|| < 1} to positive sequences in the unit ball 

Bi. 

We will show that whenever t > 0, T is a contraction on the positive sequences for which 


Cl < ||x|| < 1. Let x,y be positive sequences with ci < ||x||, |||/|| < 1, then from (2.2) 


|(Tx)o-(T2/)o| = 


(Tx)o + xi 
1 

(Tx)o + (T?/)o 


t — 


C^y)o + yi 

1 1 


+ t 


(Tx)o + Xi (T|/)o + 2/i 


From Lemma 2.3 it follows that (Tx)o > coBiit) and (T?/)o > coi?i(t), hence 

1 |(T 2 /)o - (Tx)o + 2/1 - xi| 


(Tx)o - (T 2 /)o < 


2coSi(t) [(Tx)o+ Xi][(T 2 /)o+ 2 /i] 


Lemma 2.3 also gives Xi > CiBi{t) and yi > CiBi{t) so that 

1 


|(Tx)o-(T2/)o| < 
From this we hnd 


B,{t) 


1 - 


2co(co + ciYBl{t) 

(Tx)o - (T 2 /)o 


(Tx)o - (T 2 /)o + |xi - yi\ 


For t >0, the function 


2co{co +ciyBf{t)\ Bo{t) 
1 


< 


Fi - 2/11 


2co(co + Ci)‘^Bf(t) _ 
(Tx)o - (T 2 /)o 


< 


2co{co +ci)‘^Bf{t) Bi{t) 
> 0, so that we can write 

^i(^) \xi-yi\ 


Bo{t) - Bo{t) [2ci(co + ciyBfit) - 1] Bi{t) 
The function on the right is a decreasing function on t > 0, so that 

Bi{t) ^ 5i(0) 


Bo{t) [2co(co + ci)^Bf{t) - 1] - 5o(0) [2co(co + c,)^Bl{t) - 1] 
approximately. For this reason, we have 

l(Tx)o - (T2 /)o 1 |xi - yi\ 

- < C2- 


= 0.568967, 


Bo{t) 


Blit) ’ 


C2 = 0.568967.... 


(2.9) 










































A similar calculation holds for n > 1: from (|2.1|) we find 
|(Ta;)„ - {Ty)n\ = 


n + 1 


+ 


n 


{Tx)n + Xn+i (Ta;)n + 


+ t 


n + 1 


n 


+ t 


(Tj/)n “1“ yn+1 (T?/)n “1“ 2/n—1 

n+1 n+1 


(Ta;)„ + (T?/), 


n 


{Tx)n + Xn+l {^y)n + l/n+1 

n 


{i:x)n + Xn-l {^y)n + yn-l 


Lemma 2.3 gives (Tx)„ > Cii?„(f) and (T|/)„ > Cii?„(f) so that 


|(Tx)„ - (Tj/)„| < 


(n + 1) [|(T?/)„ - (Ta;)„| + \yn+i - Xn+i\\ 

2CiBn{t) I [(Ta;)n + Xn+l] [(T|/)n + 2/n+l] 

n[\{Ty)n - (Ta;)n| + \yn-i - Xn-i\] 


[(Ta;)n + Xn-l] [(T?/)n + yn-l] 

Since Xn+i,yn+i > CiBn+i(t) and Xn-i,yn-i > Cii?„_i(t) we have (together with the lower 
bound from Lemma 2.3) 


|(Ta;)n - {Ty)n\ < 


\{Tx)r,-{Ty\ 
2clBn{t) 
x-y\\ 


n + 1 


+ 


n 


[Bn{t) + Bn+l{t)] [Bn{t) + Bn-l{t)] 


+ 


{n + l)Bn+i{t) 


+ 


nBn-i{t) 


‘^c\Bn{t) \ [Bn{t) + Bn+i{t)f‘ [Bn{t) + Bn-l{t)f‘ ) 

which can be rearranged so that 


n + 1 


+ 


n 


‘^c\Bn{t) \ [Bn{t) + Bn,+i{t)f [Bn{t) + Bn-l{t)] 


\{+x)n-{Ty)r 


Bn{t) 


< 


x-y\\ f (n + l)Bn+i(t) 


+ 


nBn-i(t) 


2cfBn(t)^ [ \_Bn(i) + Bn+lif)] [Bn{t) + Bn-lif)] 
The function on the left 


f{n,t) = 1 


n + 1 


+ 


n 


2c\Bn{t) \ \J3n{t) + Bn+l{t)\ [Bn{t) + Bn-l{t)\ 

is bounded from below by /(n,0) since Bn{t) > (2n + 1)^/^ for any n > 0. In addition, 
f{n, 0) is an increasing function of n and therefore /(n, 0) is bounded from below by its 
value at n = 1, which is approximately /(1,0) = 0.507422. Since f{n,t) > 0.507422 > 0, we 
can write 


\{+x)^-{+y\ 

Bn{t) 


+ 


nBn-i{tf 


< 


(n + l)Bn+i{t) 

[Bnif) + Bn+l{t)Y [Bnif) + Bn-lit)] 


2c\Bn{t) - 


n + 1 


+ 


n 


[Bn{t) + Bn+l{ty\ [Bn{t) + 


\x-y\\. 
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We observe that for each t > 0, the expression on the right is an increasing sequence for 
n > 1 , therefore it is bounded from above by the values as n —>■ oo, which is approximately 
C 3 = l/( 8 ci — 1) = 0.928273. Consequently, for any n > 1 


|(Tx)„-(T7/X 

Bn{t) 


<c^\\x — y\\, C 3 = 0.928273 .... 


Combining (2.9) and (2.10) then gives 


( 2 . 10 ) 


||Ti-T!/|| < 0.928273||x-i/||, 


which shows that T is a contraction. Since the unit ball with the norm || ■ || is complete, 
one can use Banach’s hxed point theorem to conclude that T has a unique hxed point b for 
which Tb = b. The sequence b = ( 6 n)n>o is positive and it is a solution of the alternative 
discrete Painleve equation (1.4). The contraction factor can be improved by improving the 


upper and the lower bounds in Lemmas 2.1 and 2.3 The lower bound in Lemma 2.3 can be 


used to get a better upper bound in Lemma |2.1[ which in turn can be used to improve the 
lower bound in Lemma 12.31 


We have not yet shown that bo{t) = — Ai'(t)/Ai(t); see Corollary 5.2 


3 Behaviour of the unique positive solution 


The unique positive solution of (1.4) necess arily satishes bn{t) > \/i, which is an immediate 
consequence of (1.4). Furthermore, Lemma 2.1 implies ^/t < bn{t) < Bn{t), for any n > 0 


and t > 0. This is illustrated in Figure 3.1 


or n = 5 and n = 10. Moreover, recalling the 




B5(t),b5(t), \/i Bio(t),bio(t), \/t 


Figure 3.1: Plots of Bn{t) (dot-dash line), bn{t) (solid line) and \/i (dashed line) for n = 5 
and n = 10 . 
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definition of B^it) given by (2.4), we observe that linit_,.oo Bn{t)/y/t = 1, and consequently, 


\im bn{t)/Vi = 1. (3.1) 

t^oo 


Remark 3.1. Observe tha t for any t > 0 and any integer n > 0, the upper bound B^V) 
for hnit) given in Lemma |2.l| satishes Vi < Bn{t) < y/f + (2n + l)^/^. For values of f > 

9 * 

3 (n + I) , Bn{t) can be written as 

Bn{t) = ^Vsicos arccos [3^'^^(n + } , 

for any n > 0. Since 0 < | arccos for f > 3 (n + it follows that 

Vi < Bn{t) < ^VVi. 

Remark 3.2. The Backlund transformations pair 


bnif) + bn+iit) 

bn{t) + bn-l{t) 


2 (n + l) 

Kit) -1 + b'Vt)' 

2n 

Kit)-t-Kit)' 


(3.2a) 


(3.2b) 


relate solutions of the equation (1.4). In fact, ( 1.4[ ) can be obtained by eliminating b'Vt) 
between (3.2a) and (3.2b). Hence, 


Kit) = 


n + 1 


n 


bnit) + bn+lit) bnit) + 6n-l(t) ’ 


which, recalling (1.1b), corresponds to 


Kit) = an+lit) - Unit). 


Moreover, (1.1a) and the latter equality yield 


bViit)-Kit) = bViit) + bVt)- 


(3.3) 


(3.4) 


The identities (3.3) and (3.4) allow us to conclude the following result. 


Lemma 3.1. For each integer n > 0, b'^{t) > 0 for all t > 0 if and only if an+i{t) > aVt), 
and this implies b^+iit) > bn{t) and also a'Vt) < 0. 


Numerically it is evident that bnit) is a strictly increasing function of t. However, it 
remains an open problem to prove this (analytically). 

Conjecture 3.2. For any f > 0, 0 < b^it) < ^/VVi)- 

Notwithstanding this, we can still prove the result for values of t 

Lemma 3.3. For each integer n > 0, we have b'^{t) > 0 for any t > 


larger than 2 

/ (2n+3)n3 \ 

y (2n+l)(n+l) y 
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Proof. Since Oq = 0, then (3.3) readily implies 6g(t) = ai(t) > 0, for any t > 0. Based on 


(3.3) together with the inequalities 

n 


< Or, < 


n 


2^t B^it) + 

which are valid for any n > 1 and f > 0, it follows that 

n + 1 n 


Kit) > 


Bfi+lit) + Bn{t) 2^/1 


Since Bnit) is an increasing function of both n and t, bounded from below by \/t and 

B (t) 

lim "I ’ = 1, it follows that for each n > 1 there exists t* such that for any t > t*, we 

t—>+oo vt 

have Bnit) + Bn+iit) < (l + ^) 2^/1, which implies > 0. 

Indeed, let t* be such that i?n+i (^n) = Bnit) is by dehnition the positive 

solution of y{y‘^ — t) = 2n + 1, it readily follows that f* must be such that 


n + 1 


[(n + 1)2 = 2n + 3, 

. Then for t > ty Bnit) + < 2Bn+iit), for each 

□ 




and hence f* = 
n > 0. 


(2n+3)n® 

(2n+l)(n+l) 


2/3 


In the same manner we can show the existence of tf such that for any t > tf, we have 

Kit) < l/i2^/t). 

Remark 3.3. The asymptotic behaviour of bn for large t in (3.1) allows us to conclude from 
(1.1b) and (3.3), respectively, that limi_,.+oo {anit) 2\/t/n) = 1 and limt^+oo {b'nit) 2\/f) = 1. 

Remark 3.4. As bnit) > y/t for any n > 0, then (3.2a) implies 

2 (n + l) 


0 < bnit) + bn+lit) = 


whereas (3.2b) implies that 


0 < bnit) + bn-iit) — 


Hence, we have 


Kit) - t + b'^{t) 

2n 

Kit) - t-b'nit)' 
> 0 , 


Kit) - t + ^bl{t) - t - b'^it) 

i.e., -iblit) -t) < b'^(t) < blit) -1. 
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4 Relation with the second Painleve equation 


Suppose we denote solutions of Pn (1.3) by w{z]a)^ then the three solutions w{z]a) and 
a ± 1 ) are related by 


a + I 


+ 


a 


w{z-, a + 1 ) + w{z; a) w{z; a) + w{z] a — 1) 


+ 2w‘^{z-, a) + z = 0, 


(4.1) 


see [m Eq. 32.7.5 on p. 730] and [3, Eq. (1.21)]. The latter is equivalent to the difference 
equation (L4|). If we set a = n + | then bn{t) = + |). It is known that 

T/^V» I P h* l—f I 4 -I \ h* O OI 1 T P 1 vn O 1 -t- O /^T A 1 V»T r T 1 T /TT- 1 O ( n Wl -(') 4 ill 11 \ 


for a = n + I then Pn (1.3) has solutions in terms of Airy functions (see [IHl 32.10(ii) on 
p. 735] or in §7.1 on p. 373]). The simplest “Airy solution” is 




0 '(^) 
0 (Z) ’ 


where 


^(z) = Cl Ai(-2-‘-'V) + C 2 Bi(-2-‘''3;), 


(4.2) 


with Ai(t) and Bi(t) the Airy functions and Ci and C 2 arbitrary constants. Observe that 
the solution w{z] \) only depends on the ratio Ci/C 2 - We now have that 

= _2-V3„(_2i/3*;l) = 2-V3^1gtA, 

and since 0(— 2 “^/^t) = Ci Ai(t) + C 2 Bi(t), then we hnd 

OiAi'(t) + C' 2 Bi'(t) 


ya{t) = 


Cl Ai(t) + 02 Bi(t) ■ 


(4.3) 


For t = 0 this gives 


6o(0) = 


Oi Ai(0) + 02 Bi(0) 
where we used the initial values 


OiA7(0) + O 2 B 7 ( 0 ) ^ 3[r(|)]" ( Ci-V^C 2 

27r Ic'i + V302 


Ai(0) = 


1 1 3V6r(|) 32/3r(i) 


32/3r(|) 


, Bi(0) = 


3V6r(|) 


, Ai'( 0 ) = 


27r 


Bi'(0) = 


271 


see [HI 9.2(ii) on p. 194]. Thus choosing O 2 = 0 gives the initial value mentioned in Theorems 
II.II and [T 2 ] for f = 0 . 

More generally, we have the solution 


/ In *^1 

w[z]n + 2 ) = q— In 


Qn{z) 


dz Qn+l{z)' 


Qn(z) = det 


di+fc 

d^ 


c^{z) 


n-l 


j,k=0 


(4.4) 


for n > 1 , with (j){z) as given by (4.2), cf. [H ITUIIT7] . 


13 




























Consider the Backlund transformations 


Vn+l — —Vn + 


Un—l Un 


2 (n + l) 

yl + y'n-i' 

2n 


yl-y'n-t' 

where ' = d/dt. Eliminating y'^ yields the recurrence relation 


(4.5a) 

(4.5b) 


n + 1 


+ 


n 


2/n+l 4“ 2/n dn 4“ 2/n—1 




(4.6) 


which is alt-dPi (1.4), whilst letting n —)■ n + 1 in (4.5b) and then substituting (4.5a) yields 

d^ 2 /r, 




= - ‘^tyn - 2n - 1, 


which is (1.2) with a = n + |, and so is equivalent to Pn (1-3) with a = n + 
Lemma 4.1. Suppose and yn satisfy the discrete system 

Xn 4“ Xyi-\-\ 2/n, 

^ni^yn 4“ 2/n—l) 


and yn satisfies (4.5). Then Xn and yn satisfy the system 

dx„ 

= -2xnyn + n, 


dt 

dyn 

dt 


dn t, 


and Xn satisfies 


d'^Xr, 

dt^ 


2 x,. V^> + 


n 

2Xr. 


Proof. From (4.8b) and (4.5b) we have 

2n 


n 2 dyn 

= 2Xn = yn -^ - C 


dn 4“ 2/n—1 dt 

from which we obtain equation ( |4.9b ). Then differentiating (4.9b) gives 

d^yn „ dyn AXn . 

—^ = 2yn— - 2^- - 1 . 


dt" 


dt 


dt 


(4.7) 


(4.8a) 

(4.8b) 


(4.9a) 

(4.9b) 

(4.10) 


Substituting for the derivatives of yn using (4.9b) and (4.7) yields equation (4.9a). Finally 
solving (4.9a) for yn and substituting in (4.9b) shows that Xn satisfies equation (4.10), as 
required. □ 


14 














































We remark that making the transformation = —2 with z = 


eqnation (4.10) yields 


dz^ 2v \dzJ 


zv 


n 


which is known as P 34 , as it’s eqnivalent to equation XXXIV of Chapter 14 in 


The “Airy-type” solutions of (4.7) and (4.10) respectively have the form 


a\ '^n{t]0) 

yn{t; 0) = — \n 


dt r„+i(f; 6 ')’ 


d^ 


in 

(4.11) 


(4.12) 


where 


= det 


di+fc 

df^ 


n n—1 




(4.13) 


j,k=0 


with To{t-,6) = 1 and (fit] 6) = cos( 6 *)Ai(f) -|-sin( 6 *) Bi(f), with Ai(f) and Bi(t) the Airy 
functions and 9 G [ 0 , tt) a parameter; we have set Ci = cos( 6 *) and C 2 = sin( 6 *) to reflect that 
the solution only depends on the ratio of the constants. The “Airy-type” solutions (4.12) 
are derived from the “Airy-type” solutions of Pn (1.3) given in [U UHl ITTj . 

The simplest non-trivial “Airy-type” solutions of (4.7) and (4.10) respectively have the 
form 


yo{t; 9) = 
xi(f; 9) = 


■—In 9) = - 
~\n^{t-,9) = 


cos{9) Ai'(f) -|- sin( 6 ') Bi'(f) 
cos( 6 *) Ai(f) -|- sin( 6 *) Bi(f) 
cos( 6 ') Ai'(f) -|- sin( 6 ') Bi'(f) 
cos( 6 ') Ai(f) -|- sin( 6 ') Bi(f) 


t, 


(4.14a) 

(4.14b) 


recall that Xo(t; 9) = 0. The structure of the solutions (4.14a) and ( |4.14b ) depends critically 
on whether the parameter 9 is zero or not which is shown in Lemmas 4^ and m below. As 
Ai(t) decays exponentially as f —)■ 00 , whereas Bi(t) increases exponentially as f —)■ 00 , so if 
9^0 then Bi(t) will dominate for large positive t. 

If we seek a solution of (4.7) with ynit) as f —)■ 00 , then necessarily c = ± 1 . 

Hence the following asymptotic series are easily derived 


y, 


-{t) 


2n + 1 I2n^ + 12n + 5 


At 


1/9 -|- 1 

ynit) =-t"'" + + 


32f5/2 

I2n^ + 12n + 5 


32t5/2 




(4.15a) 

(4.15b) 


It should be noted that there are no arbitrary constants in these asymptotic expansions, 
they occur in exponentially small terms. 


Lemma 4.2. IfyQ{t]9) is given by (4.14a), then 


yoit;9) 


+0(^-1), 
-t^/^ + o{t-^), 


if 9 = 0, 
if 9y9 0. 


(4.16) 
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Proof. Using the known asymptotics of Ai(t) and Bi(t) 


Ai(t) 

Ai'(t) 

Bi(t) 

Bi'(t) 


.-C 


UAe-C 


20 F 


1 + 


fiAeC 


1 + 


48 13/2 

7 

48 PP 
5 


+ 


+ 


TT 


1 - 


48t3/2 

7 


+ 


385 
4608 P 
455 
4608 P 
385 


85085 




+ 


663552 

95095 

66355219/2 ^ ^ 


4608 

455 


+ 


85085 


66355219/2 
95095 


+ o{t-^) 


48 ^3/2 4608 P 663552 ^9/2 


+ o(r®) 


with ( = see [IHl 9.7(ii) on p. 198], then if 6 * = 0 , as t —)■ oo 


!/„((; 0 ) = -AU = (V 2 


1 + 


7 


455 

+ + 


95095 


48f3/2 4608 f3 663552 f9/2 


+ o{t-^) 


Ai{t) 


1 - 


5 385 

+ 


85085 . _e. 

48f3/2 ' 4608 f3 663552 f9/2 ^ ^ 


1 


= tV 2 + A 


+ 


15 


4f 32f5/2 64 


+ 0(riV2) 


whilst if 0 ^ 0 , as f —)■ cxD 

cos( 6 ') Ai'(f) + sin( 6 ') Bi'(t) 


yo{t-,0) = 


cos( 6 *) Ai(f) + sin( 6 *) Bi(f) 
= -tV 2 


1 cos(6') e j 

4 +48(3/2+C’(‘ )j 

|- + sin(6') e'^ j 

48f3/2+^(^ )} 

1 cos(6') e”*’ j 

48f3/2 )j 

|- + sin(6') e’’ j 

d + ' ‘2(‘ )} 


= -tV 2 


1 - 


7 


455 


95095 


48^3/2 4608 f3 663552 f9/2 


+ o{t-^) 


1 + 


48 f3/2 


+ 


385 85085 . .g. 

4608 f3 ^ 663552 ^ ^ ^ 


+ 0(fB2e-2C) 


= _tV 2 + ^ ^ 

^4f^32t5/2 


+ 


15 

64 P 




□ 


Plots of ?/o(U and Xi{t; 6) for various values of the parameter 6 are given in Figure 4.1 
We note that 


2/o(O;0) 


cos( 6 *) Ai'(O) + sin( 6 *) Bi'(O) 
cos( 6 *) Ai(0) + sin( 6 *) Bi(0) 


33/6 A - Vs tan 0 \ . 2 
^ l^l + VStan^; ^ 


which is a decreasing function for 0 < 0 < A- 
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Figure 4.1: (i) Plots of yo{t;9) for 6 ^ = 0 (upper curve), 9 = Itt, ^vr (lower 

curves, with 9 = the highest and 9 = the lowest), and the curves y = ±\/t (dashed 
lines).(ii) Plots of Xi{t]9) for 0 = 0 (upper curve), 9 = ^vr, ^vr, ^tt, ^tt, (lower curves, 
with 9 = ^TT the lowest and 9 = ^tt the highest), and the curves x = ±l/{2\/t) (dashed 
lines). 


An analogous situation arises for the solution yn{t]9) given in (4.12), i.e. 


= P + .f » = o. 


(see Lemma 5.1 below). We remark that Fornberg and Weideman P, Fig. 3 on p. 991] 
plot the locations of the poles for the solution w{z-^ |) of Pn (1.3), which is equivalent to 
1/2 (t; 0), for various choices of 9. These plots show that the pole structure of the solutions is 
signihcantly different in the case when 0 = 0 compared to the case when 0 7 ^ 0 . A study of 
“Airy-type” solutions of Pn ( |1.3[ ) and P 34 (4.11), which are equivalent to //n(t; 0) and Xn{t; 0), 
for various choices of 0 is given in j5]. 

Lemma 4.3. Ifxi{t;9) is given by ( 4.14b[ ), then 

xi{t]9) = J 2 

Proof. The proof is very similar to that for Lemma |4.2| so is left to the reader. 


1 ^ 1/2 + 0 (^- 2 ), 
. 1 ^ 1/2 + 0 (^- 2 )^ 


if 0 = 0 , 
if 0 ^ 0 . 


(4.17) 


□ 


17 




















5 Airy solutions of alternative discrete Painleve I 


an(t) =Xn(t-,0) = bn(t) = 0) = — In 


Now we consider the case when 6 = 0, i.e. 

A' 

dt 

where 

An(t) = rn{t;0) = det 


A„(t) 


dt A„+i(t)’ 


Aj+k 1 n 1 


(5.1) 


(5.2) 


with Ao(t) = 1. We remark that A„(t) given by (5.2) arises in random matrix theory, in 


connection with the Gaussian Unitary Ensemble (GUE) in the soft-edge scaling limit, see 
e.g. [ini p. 393]. Further, for n > 1, A„(t) has the multiple integral representation 


A„(f) = 


(- 1 )' 

(27ri)’ 


exp - txj) {xk - xif dxi ... dxn, 




l<k<l<n 

in which C is the standard Airy contour from cx)e“^'A ^ cxie^'A p. 393]. 


Lemma 5.1. If an{t) and bn{t) satisfies the recurrence relation (1.1) with 

ao(t) = 0, bo{t) = -Ai'(f)/Ai(f), 
where Ai(t) is the Airy function, then as t ^ oo 

for n>l, 

bn{t) = t^^^ + 0(t~^) , for n>0. 


(5.3) 

(5.4a) 

(5.4b) 


(4.5a). 


Proof. We shall hrst prove (5.4b) by induction. If bn{t) = yn{t', 0) satishes the recurrence re¬ 


lation (4.6), which is a consequence of (1.1), then 6„(t) also satishes the differential-difference 


system (4.5) and the differential equation (4.7). We shall use the Backlund transformation 


Glearly (5.4b) is satished for n = 0 from Lemma 4.2 Now we assume as the inductive 


hypothesis bn{t) = 0(t ^), so we know from ( |4.15a[ ) that 

^ ^ 1/2 , 2n + l 12n2 + 12n + 5 , (2n + l)(16n2 + 16n + 15) , 

On{l)-l + qo+ 5/2 + fi/l+4 ), 


dt 32^5/2 

and hence it is easily shown that 


, 2 ,, 2n + l 2n‘^ + 2n + l 5(2n + + dn + 5) 

^^_ 2 ^ 5(12n^+12^ 

2U/2 4t2 + 64U/2 )■ 


+ 0{t-^) 


blit) + blit) - t = 


n -|- 1 
tV 2 


n -I- 1 5(4n^ + 8n + 5) 


2f3/2 


+ 


32 
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Therefore 












































then 


2(n + l) 


blit) + b'^it) - t 




n + 1 12n^ + 24:12 + 17 


and so from (4.5a) 


bn+iit) — —bnit) + 


t 


2(n + l) 


16t5/2 


+ o{t-^), 


blit) + h'^it) - t 


= + 


2?7, + 3 12n^ + 36?7, + 29 


4t 


32f5/2 




which is (4.15a) with n —)■ n + 1. Hence the result (5.4b) follows by induction. The result 
(5.4a) then follows immediately from (1.1b). □ 


Corollary 5.2. The positive solution bnit) in Theorems 1.1 and 1.2 corresponds to the initial 
value boit) = — Ai'(f)/Ai(f). 




(ii) 


Figure 5.1: (i) Plots of bnit) = Unit', 0) for n = 0,1, • • •, 6, and the curve y = \/i (dashed 
line), (ii) Plots of anit)/n = Xnit]t))/n for u = 1,2, ..., 6 , and the the curve x = l/i2\/t) 
(dashed line). 


Plots of bnit) = ynit;0), for u = 0,1,..., 6, and plots of a„(f)/n = a;„(f;0)/n, for 
n = 1,2,. .., 6 , are given in Figure |5.1[ These plots suggest the following conjecture, which 


corroborate Conjecture |3.2| and Lemma 3.3 
Conjecture 5.3. If 0 < ti < t 2 then 


bnitl) < bnit2), ttnih) > a„(t2), 


(5.5) 


i.e. bnit) is monotonically increasing and a„(f) is monotonically decreasing for t > 0. For 
fixed t with t > 0 then 


Vi < bnit) < bn+iit), 


1 ^ ttnit) ^ ttn+lit) 

2Vi n n + l 


> 0 . 


(5.6a) 
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